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Abstract 

We propose a method for extending a given asset pricing formula to account 
for two additional sources of risk: the risk associated with future changes in 
market-calibrated parameters and the remaining risk associated with idiosyn- 
cratic variations in the individual assets described by the formula. The paper 
makes simple and natural assumptions for how these risks behave. These extra 
risks should always be included when using the formula as a basis for portfolio 
management. We investigate an idealized typical portfolio problem, and argue 
that a rational and workable trading strategy can be based on minimizing the 
quadratic risk over the time intervals between trades. The example of the vari- 
ance gamma pricing formula for equity derivatives is explored, and the method 
is seen to yield tractable decision strategies in this case. 



1 Introduction 

Any rational and practical approach to managing a financial portfolio naturally involves 
finding an optimal balance between a number of competing criteria. For example, a 
manager will wish to make good (but not excessive) use of available market data. They 
will also need to work with a pricing model in order to estimate risk and project to 
the future. However, the pricing model must balance complexity and accuracy against 
simplicity and implementability. They will want to trade but not too frequently, and 
before each trade will need to choose from many alternatives. They may need a set of 
investment decision strategies which are explicable to their boss, and will stand up to 
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stringent scrutiny. They must be aware that strategies and parameters ahke need to 
be updated continuously. 

The paper focuses on an ideahzed market of derivatives on a single underlying asset 
St over a time period t e [0, T]. We address the hypothetical situation of a manager 
who has written (sold) to a client a large and risky over-the-counter (nontradable) 
contingent claim F on S, with maturity date T (which we think of as quite large, 
say three years). Having accepted the large fee, and the large liability, the manager 
must now create a portfolio of investments in the exchange traded derivatives with the 
aim of adequately or even optimally hedging the risk. Thus the market consists of a 
number of different securities: a bank account, a non-dividend-paying stock, a family 
of derivatives on the stock and finally the over-the-counter claim. We assume 

1. the bank account pays and charges a constant rate of interest r, which we take 
to be zero for simplicity. There is no limit on borrowing; 

2. the stock price at time t, denoted St, is always positive; 

3. the exchange traded derivatives are taken to be european puts and calls over 
a variety of strikes K"' and maturities T° whose prices at time t are denoted 

a — 2, . . . , M. By convention, we denote the stock itself by D"^ with a — 1. 

4. the over-the-counter claim price ft is the value at time t of a claim given by 
a function F of the stock price history {St}t<T- That is. Ft is an J^t random 
variable. By convention, we denote the over-the-counter claim value by D'^ with 
a = 0. 

5. the market is frictionless: there is perfect liquidity, no transaction costs, unlimited 
shortselling, and the market is open at all times. 

We assume the manager will set up an initial portfolio att — O, and then be making 
single trades at certain times ti < t2 < ■ ■ ■ < tk < ■ ■ ■ < T. In principle, these will 
be stopping times, but for this paper we take the deterministic values ti = iSt for St 
fixed. Our problem is to provide this manager with a consistent, workable strategy for 
trading, and to measure the performance of this strategy. The following two questions 
are addressed in this paper: 

Ql What is the condition for optimal allocation at time t — 01 

Q2 What is the optimal single (or double or triple) trade which can be made at time 
ti,i > 0? 

It is supposed that the rational manager will base their strategy on a derivative 
pricing model whose historical performance has been stiidicd. bcnchmarked, and found 
to lead to acceptable modeling errors. This model is taken to be a finite parameter 
"risk neutral" pricing formula, whose parameters 6*", a = 1, . . . , are to be calibrated 
to the observed prices of the M — 1 derivatives (where N « M). Of course, in 
contradiction to the modeling assumptions, the calibrated values 9^ will depend on the 
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time t, reflecting changing market conditions. Furthermore, the cahbrated model will 
generally fail to match all observed prices at a given time, the differences being thought 
of as "idiosyncratic errors" . The key point of the paper is to extend the original model 
by including minimal additional assumptions on changing parameters and idiosyncratic 
errors which are consistent with the observed performance of the model. These extra 
effects must be included in any rational portfolio strategy. 

We argue that a rational strategy for the manager will be to minimize the total risk- 
neutral variance of the portfolio returns over each inter-trading time interval, without 
regard to the mean portfolio return. "Risk-neutral" refers to the measure which cor- 
rectly prices the current market but not historical values. Variance optimization is 
equivalent to optimization with a quadratic utility function, which is an acceptable 
approximation to the general utility function over short time-intervals. Disregarding 
the mean return is consistent with a no-arbitrage condition true in an efficient market, 
but also justifiable pragmatically because the manager is profiting from the large fee 
collected for underwriting the claim F, and needs only consolidate that fee with mini- 
mal risk. Furthermore, estimating mean returns over short time intervals is a game for 
speculators and arbitrageurs, not hedgers. 

Before each trade the manager will perform a number of steps: 

1. place the payoff of newly expired options (if any) in the bank account; 

2. observe the market prices; 

3. recalibrate the risk-neutral parameters to the new data; 

4. estimate the value of the liability F using the recalibrated pricing model; 

5. find the optimal single trade which will minimize the estimated risk-neutral vari- 
ance of the portfolio over the period until the next trading time. 

The main theorem of this paper states that the risk-minimal portfolio for any trading 
time-interval is unique, as is the optimal single trade at any time. 

We then discuss models such as the VG model ||MCC99|| for which the stock process 



is Markovian and has a closed form characteristic function. It is shown that the Fourier 
transform pricing method of ||CMOCI|| extends and allows for efficient computation of 
portfolio variances, and hence optimal portfolios and trades. 

One important outcome of the present discussion is a clear operational definition 
of essential concepts of hedging which shows how these concepts can be refined and 
extended: delta hedging, gamma hedging, vega hedging, model recalibration. Further- 
more, the inclusion of risk associated to "idiosyncratic errors" breaks the degeneracy 
(nonuniqueness) of the hedging decisions derived from naive delta-hedging (in practise 
by biasing trading toward near-the-money instruments). For a detailed discussion of 
hedging and derivatives from the trading perspective, see ||Tal97|| . For a discussion of 



model risk, see [ HebOl ]. For another approach to optimal investment in derivatives see 
| CJM01|| . 



4 



2 The model 

2.1 The option pricing formula 

The option pricing formula is assumed to arise via arbitrage pricing theory set in the 
risk-neutral filtered probability space {Q,J-',{J-'t}te[o,T],Q.) over a time horizon [0,T]. 
For each fixed set of parameters 6 = {O"-)^^^ G R^, the stock is a positive Q-martingale 
("cadlag") (recall the interest rate is taken to be zero). We consider models in which 
St is one component of a multidimensional Markov process {St, Yt) whose values are 
observable in the market: in our examples Y will have dimensions zero and one. Let 
yt = (T{Sr, Yt-, t < t} G J-'t he the filtration of market observables. 

The pricing formula giving the value at time t for any European style contingent 
claim with payoff Ft at date T is simply 



F{t, S, Y, 9) = EiFTiSrm, St = S,Yt = Y,9 fixed), t e [0, T] 



(1) 



Examples: 

1. Black-Scholes model: The single parameter can be taken as ^ = logo" and the 
stock process is given by St = Soexp[e^Wt — (there are no extra Y vari- 

ables). The pricing formula if Ft = {St — K)^ is simply the Black-Scholes call 
option formula. 



2. Variance-Gamma model ( ||MCC99| ) Here the log return process Xt = log (St/So) 
is defined to be the Levy process | |Ber9(j[ |,[ [lW89|| : 



Xt = fi0t+ /* n xN^'''>\dxdt) 

Jo J-oo 



(2) 



whose jump intensity measure is 



\x\/r]± 



ueix 



\x\ 



±x > 



(3) 



for three positive parameters 6 = {a,ri^,ri_). The characteristic function turns 
out to be 



<l>xAu) = E(e 



iuXt 



(4) 



The martingale condition on S needs ?7+ < 1 and is then equivalent to $(— «) = 1, 
so 
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With this condition satisfied, then the stock process St has the form 



St = So + 



l)Sr-N^'''\dx dr) 



where N denotes the compensated (martingale) process 

iV(^«)(da; dt) = N^'''\dx dt) - ueix) dx dt 



(6) 



(7) 



3. Stochastic volatihty models: Standard stochastic volatility models |Pes93|] ,||PS001 
|PS00a|| , piF50q] such as 

V^t St dWl 



dSt 
dvt 



(a - bvt) dt + a ^Jvl dWi 



(8) 



can be considered. This model is Markovian in {St,Vt), not in St alone and the 
option pricing formula at time t depends on the values {St, Vt). Our approach will 
be to follow |[BZ01|| and treat the stochastic squared volatility Vt as the "effective 
observable" defined by 



Vt 

where Si = St-{N-i)At/N- 
d{W\W^) = pdt. 



N 

lim lim N-^ y\\og(Si/Si^i] 



(9) 



This model has four parameters (a, b, a, p) where 



Many models like the VG model and some stochastic volatility models have a closed 
formula for the characteristic function In such cases ||CMOO|| have shown how the 
Fast Fourier transform method provides a numerically efficient method for evaluating 
(0). For example, the European call payoff function with a strike K = can be written 
as a complex contour integral along the shifted contour {—oo — ie, oo — ie) for any e > 0: 



where 



1 

27r J-c 
Q{u) 



S 



(10) 



'111 



tu — 

To take the expectation to evaluate (0), plug in this formula and interchange the du 
and E{-) integrals 



F(t, S) 



1 rco—ie . 

— / E e^("- 



2tx 
1 

2^ 



s. 



iu+1 



^{u - i) 



T-t 



™^Q(m) du 



du 



'12) 



The explicit integral can now be evaluated (approximately) very efficiently for a linearly 
spaced family of log strike values {ko + n5k}n=Q,±i,±2,... using the fast Fourier transform. 
Similar formulas hold for put options and the over-the-counter claim itself if its payoff 
has a known Fourier transform. 
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2.2 Varying market conditions 

A careful trader using a pricing formula such as (|1|) will wish to correct for changing 
market conditions. It will be supposed that prior to its practical application (to hedg- 
ing) the performance of the formula has been carefully benchmarked against real market 
data and a historical time series of mean square estimated values 9t (as described in 
section §3.1) has been observed. When we allow the vector of parameters 9 to become 
a process, it is natural to pick a gaussian mean-reverting process and therefore we 
assume: 

Al The time varying parameters 91 form an R^-valued Ornstein-Uhlenbeck Ito dif- 
fusion (gaussian, mean reverting). The noise process driving 9t is independent of 
the noise process driving St- 

Let = a{9s,s< t} C J'f 

A consistent pricing formula based on the filtration yt x is given by 

F{t, S, Y, 9) = E{FT{ST)\yt X Gu St = S,Yt = Y,9t = 9), t G [0, T] (13) 

Note that F depends on the detailed modeling of the process (described of course by 
extra parameters), and hence F ^ F. For the present paper we regard use of F in place 
of F as a change in the underlying pricing model, and is hence "against the rules" . 

Since our purpose is rather to use F to find a hedging strategy which will include 
low order approximations to account for variations in 9, the point of view of varying 
parameters can be simplified. We adopt a natural compromise position which is to 
regard the parameters 9 as constant over the intertrading intervals [)f:j_i,tj). At each 
time ti the parameters are assumed to jump according to a discrete R^-valued OU 
process which is independent of the stock price St^. Over the time period [tj_i,tj), 
the stock price is assumed to evolve according to the model with fixed 9 = 9t^_^. The 
distribution of 6'^. conditioned on 6'(._j is gaussian. To keep the discussion simple, we 
will focus on the conditional covariance estimated by the quadratic variation statistic: 

1 " . 

0"' = -E[^|-^LJ[^i-CJ' -.h=l,...N (14) 
fc=i 

obtained from a series of times tQ^ii, . . . ,in (all prior to the trading times) with tj — 
ti_i = St. 

2.3 Idiosyncratic variations 

Even if the model parameters 9 are confidently known at time t, the observed option 
prices D°' will not fit the formula perfectly. There are certainly model errors which are 
idiosyncratic to each exchange traded security, and a careful trader must also account 
for this. As before, it is supposed that these idiosyncratic errors have been benchmarked 
by observing how the formula matches real data and that the process 

£^ = Df-F''{t,St,9t) (15) 

has been sampled over time. We assume 
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A2 the errors St form a multidimensional Ornstein-Uhlenbeck Ito diffusion, indepen- 
dent of the noise processes driving St and 9t- For simplicity, we assume that the 
covariance of Et^ conditioned on £t^_i is diagonal. 

Trading instincts suggest that the predicted variances will depend most strongly on 
two variables, namely moneyness k = K/ S and time to maturity t = T — t, both of 
which vary in time for a specific derivative a. Therefore we adopt the following crude 
method for predicting the idiosyncratic variances: Let the variance of ^^^+5^ conditioned 
on £^ be 

Vt" = S^viK'^/St, T" - t) (16) 

where is the following historically observed quantity which depends on moneyness 
and time to maturity: 

^(-^ -) = ^tT. s^is^^i^r' - ^t'^'^'r (17) 

where {tfc} is a sequence of times prior to the trading times U. Here k, r) is defined 
to be that derivative whose characteristics (-ft'", T") at time t most nearly match (k, r), 
i.e. the a which minimizes 

d{a- (/€, r)) = {K^/St - nf + ((T" -t)/r- if (18) 



3 Optimal trading 

We now develop in detail the problem of optimal trading at time tj. In what follows 
we use the shorthand notation Si = St^, etc. for quantities evaluated at time U. The 
formulation we will adopt will take into account the following information at time ti. the 
observed market prices , the empirical covariance G'^'', and the empirical variances 

We assume that immediately following the previous trading time the (self- 
financing) portfolio consists of t^'^_i units of each derivative, the liability for the claim 
Z)° (so 7r° = —1 always), and the remaining value -Bj-i in the bank account. The 
portfolio value X is 

M 

X,_i = Y.^tiDti + B;-i (19) 

a=0 

= X,„i-5:<iDti (20) 



3.1 Calibration to observed prices 



At time tj, we compare the observed market prices Df,a = 1,...,N to the model 
prices to update the market parameters 6. Following [ [MCC99|] , the vector of estimated 
parameters 6i is the minimizer of 



min{^ \\ogiD^)-\ogiF^iU,S.,e))\'} 



(21) 
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Note that by definition = and = = S so the sum is over a = 2, . . . , A'". 
Refining (pT]) by weighting by trading volume might be considered. 

3.2 Marking to market 

Given the recahbrated market parameters di and current market prices D°',a > 1, the 
updated portfoho value can be calculated for time tj. The best estimate of the current 
value of the contingent claim is the expected value at ti of the (discounted) final 
claim: 

D^ = F\U,S,J,) (22) 

(this explains why = always). Therefore the updated portfolio value, just prior to 
trading, can be written 

X, = Y: vrr_iF"(t„ S., e,) + Y: <i^r + 5.-1 (23) 

3.3 Finding the optimal trade 

To determine the optimal trade (i.e. to find the best new portfolio weights vr", Bi), we 
project 

Xi = Y: <F'^iU, S,, 9,) + J2 + Bi (24) 

a>0 a>2 

forward to the next trading time tj+i and consider the (stochastic) value Xj+i. This is 
characterized by the values tt", B^, the new market prices -D^i, and the new parameter 
values 6'j+i. Consider the difference AX = Xj+i — Xj = J2a>o tt" AX° where 

AX" = [F'^{ti^i, Si+i,9i) — F°'{ti, Si,9i)] 

+ 'S'i+i, 6'j+i) — F"(tj+i, Si+i, 9i)] + [£°'_^-^ — £°'] 

= AX^ + AX^ + AX° 

The middle term AX2 will become a bit awkward, so we reorganize it by Taylor 
expanding in powers of A9i = 9i^i — 9i: 

F°(tj+i, S'j+i, ^j+i) — ^"(tj+i, Sj+i, ^j) = ^(9g F°)(tj+i, 5'j+i, ^i)A6'i (25) 

n>l 

We have adopted a condensed multiindex notation so that = dgl . . . O^n and A9 = 
{9-^) \ . . (9^) ^ with n = (ni, . . . , n^r). n > 1 is shorthand for I]i > 1- 

In what follows, let Ei{-) denote conditional expectations EQ{-\J-'t-). We consider 
first the mean of AX. Note -E'j(AX") = since F" is a martingale. In general, 
however, i?j(AX^) 7^ 0,-E.j(AX^) 7^ since conditional means of OU processes are 
never identically zero. Statistically, mean returns are small, difficult to measure and 
very unstable in time. From a trading point-of-view, seeking mean returns over short 
time intervals is the job of the speculator, the day-trader and the arbitrageur. For 
our hedger, it makes no sense to trade on expected returns. Therefore, from both the 
statistical and trading points of view it makes sense to introduce an extra assumption 
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A3 The trader will ignore the mean returns of AX2, AX3, and will trade as if (AX^) = 
and Ei{AX^) = 0. 

Thus we assume that when conditioned on J^i, A9i and S°^^ — £^ are taken to be 
mean zero gaussians, with conditional covariances 0"'' and V°'. 
Under these assumptions, the variance or total quadratic risk 

R = veiTi{AX) = Yl K^iR"'^ 

decomposes into three terms 

R'^ = Covi(AXf; AXf) + Cov,(AX2"; AX2^) 
+ Cov,(AX3-;AXf) 

= Rf + Rf + Rf 

• Ri is the risk associated with changes in the underlier S; 

• i?2 is the risk associated with changes in the market reflected in evolving param- 
eter values; 

• R3 is the risk associated with deviations between observed prices and model prices. 

This type of decomposition is quite general. 

We have now succeeded in expressing the quadratic portfolio risk in terms of the 
pricing formulas F°, the conditional covariance matrix ©"^ and the variances V°'. 

Proposition 1 Consider a pricing model of the above type under assumptions A1,A2,A3. 
Then 

1. there is a unique portfolio tt* with tt? = —1 which minimizes the quadratic risk 

K = arg min,^,o=_i E 7r"7r^[i?f + Rf + Rf] (26) 

a,/3>0 

over the period given by the solution of 

f>oa ^ ^ R^f^T^f^ a>0 (27) 

/3>0 

2. for any portfolio Hi-i with ti^ = —1, there is a unique optimal single trade con- 
sisting of X* units of the ^th derivative, 7* > 1 where: 

(A*, 7*) = arg min,,,>, ^ + XS^) {nl, + A<5?) [Rf + Rf + Rf] (28) 

a,/3>0 

Here arg min denotes the solution of the associated minimization problem. 
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Proof: We need only note that Ri is positive definite while i?2, R3 are positive semi def- 
inite matrices. 

□ 

Selecting the optimal single trade given 7rj_i is a simple matter given the matrix R. 
The optimal amount A*'''' to trade of a given derivative 7 is 

A*'^ = (29) 



Such a trade decreases the unimproved portfolio risk J2a,i3 by the amount 
{J2i3 R^^ T^i-iY / R^'^ the optimal asset to trade is simply that 7* which shows the max- 
imal improvement. 

Similar arguments give one the optimal double, triple trade etc. 



4 Implementing the method 

Implementation of the method for a particular trading formula depends on efficient 
determination of the risk matrix i?"^ before every trade. We demonstrate how this can 
be done by generalizing the Fast Fourier Transform (FFT) method of ||CMOO|| when 
the log-return process Xf = log(5't/5'o) has a known characteristic function We 
suppose that for each a = 0, 1, . . . , M the payoff function F'^ISt) has a Fourier formula 
similar to (|T0|) but with the function Q{u) replaced by Q°'{u). 
First note that R3 is always simply the diagonal matrix 



R 



a/3 



ra/3 ci2 



(30) 



Ri, R2 are more complicated. We note that they are symmetric matrices which fall 
naturally into block submatrices 



R 



Roo 


Roi 


Roc 


Rop 


Rio 


Rn 


Ric 


Rip 


RcO 


Rci 


Rcc 


Rep 


. Rpo 


Rpi 


Rpc 


Rpp _ 



(31) 



where Roc is 1 x Mc,Rcp is Mr, x M„ etc. where M^, M„ are the number of traded calls 



and puts. Put-call parity leads to a number of relations amongst the components: 



Rrr — R. 



■cp 



R 



Rep Rpp 



R 



Ipl 



etc. 



(32) 



We indicate here how to compute the parts of Rcc with fixed maturities T", and any 
log strikes k'^^k^ . 

For Ri terms note that 



,a/3 
cc 



E, {Ft,,F^ - F^Fi^ 



(33) 
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The usual approach is to Taylor expand in powers of At which leads to generalized 
delta-gamma hedging. However with the FFT method, we can perform the one re- 
maining expectation leading to an explicit double Fourier integral 

Da/3 _ 1 
^l,cc — 



(27r) 



Q{ui)Q{u2) 



(34) 



A single application of the two dimensional FFT solves this problem for a family of 
calls with linearly spaced log strikes but with fixed dates T",T^. 
Now we seek a similar formula for 



r>al3 _ 771 
^2,cc — 



(35) 



When we plug in the Fourier integral, we highlight the ^^-dependence of $ by writing 
$(m;6'). Then using the independence of S and 9 expectations leads to 



xE, Stt^""'^^^ 



X 



^U2-i;9i+,f'-''^' -^U2-i;9 



\dui du2 



{27iyJ J' 



-iuik°' —iu2kl^ 



Q{ui)Q{u2) 



xE, [$(«i - I- 4+1)^"-*'+^ - *(«i - ^; K 



T°'-U. 



xeAsii\-^' 



^U2-i;9i+i) 



$(m2 - Oi) 



dui du2 (36) 



which becomes 
1 



s: 



iiui+U2)+2 iuik'^-iU2kf 



Q{ui)Q{u2) 



$(mi + U2- 2i; ^,)^*E, [[$(mi - i; 9,+^f - ^u, - i; 9^^-'^+' 

[^U2 - i; 9i^if'-''^' - Hu2 - i; ^0^'"*'+^]" 
^u, - %■ ^,)^*$(«2 - ^; 9,)^'Ei [$(«i - %■ 9,+if - $(mi - I- 4)^"-*'^ 



X 



xE, 



^{U2 - i; 9i^if'-''+' - ^{U2 - i] 9if'-''+''^ dui du2 



(37) 
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The remaining ^^-expectations will be quite complicated since they involve differencing 
the function A pragmatic approach is simply to Taylor expand in powers of A^. 
Then the leading term is 



AT 



a,b=l 



with higher order terms given by more complicated integrals. These integrals are effi- 
ciently calculated by FFT as before. 



4.1 Example: the Black— Scholes model 

The Fourier technique above applies to this model, since the characteristic function is 
simply 

= exp[-e2^(M^ + iu)/2] (39) 

The matrices Ri,R2,R3 can be interpreted in the standard language of hedging. 
For example, by taking the second order Taylor approximation to Ri in powers of At 
and setting Y.a,p Ri^'^°''^^ — 0, we obtain the conditions for a delta-gamma-hedged 
portfolio. Note that in our present formulation, the optimal solutions make this term 
small but not zero. Similarly, setting the lowest order At term of i?2 equal to zero 
yields a vega-hedged portfolio which shows vanishing first-order sensitivity to changes 
in the estimated volatility. Again our optimal portfolios will be approximately but not 
perfectly vega-hedged. Finally, making R^ small amounts to choosing preferentially 
amongst the more liquid derivatives, since their idiosyncratic risks are generally rela- 
tively smallest. In practise, this means favoring investment in near-the-money assets 
to far in/ out-of-t he-money derivatives. 

We illustrate the method with an example which shows the optimal hedge for a 
Euro-call contract with strike 1 and maturity T = 3 (all times in years) whose total 
Black-Scholes value is $1. In a pure Black-Scholes market where the current stock price 
is $1, with puts and calls with maturity 1/2 and 11 log strikes k = (—0.5, —0.6, . . . , 1.5), 
we calculated the optimal hedge (excluding the underlier S) for the period dt = 1/12. 
For illustrative purposes, we take the other parameters to be = 10~'^,Vc = Vp = 
10~^(6, 5, 4, 3, 2, 1, 2, 3, 4, 5, 6). The graph shows the value bought of each call option (+) 
and the value (x) sold of each put. The ratio of quadratic hedged risk to unhedged risk 
is 0.0319, the hedged to unhedged delta ratio is —0.0441 and the hedged to unhedged 
gamma ratio is —0.7411. With optimized choices for FFT parameters the computation 
required 21.3 megaflops. 
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optimal portfolio to hedge long term eurocall 



~i 1 1 1 1 r 



~i 1 1 r 



0.16 



0.14 



0.12 



? 0.1 
o 



g'O.OSf 
o 

OJ 

5 0.06 [ 



0.04 



0.02 



+: value of calls bought 
x: value of puts sold 



-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 

log strike 



4.2 Example: the VG model 

We begin by noting that for this model the parameterization 6 = {a,r]^,ri_) is not a 
hnear space, but can be associated with a natural riemannian manifold with metric 
given by the Fisher information metric. Therefore, we replace assumption Al by the 
natural family of Ornstein-Uhlenbeck processes which constitute a preferred class of 
mean-reverting processes on a general riemannian manifold. However, we expect this 
nontrivial geometry to affect only the higher order terms in the Taylor expansion for 
i?2 and for a preliminary implementation we would work in the tangent space and keep 
only first order Taylor terms. Thus we would use (^) with the estimated B matrix as 
before. 

Now, since the characteristic function for the VG model has the simple analytic 
expression (^, the important double integrals can be efficiently calculated via the 
FFT. 

It is now known that the VG model does not work well in pricing derivatives across 
maturities, and so our framework cannot be expected to perform adequately. Work is 
in progress to extend the VG model to include nontrivial time correlations. 
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4.3 Example: the affine stochastic volatiUty model 

In the model given by equation (H) the characteristic function m) = i?(e*^*") for 
Xt = \og{St/So) depends on the current level of squared volatility Vq: 



where 



u) = exp[C{t, u) + D{t, u)vo] 



C{t,u) 
D{t,u) 



{h — ipau + d)t — 2 log 



b — ipau + d 



-,dt 



dt 



g 
d 



b — ipau + d 
b — ipau — d 

\J {ipau — bY — a'^{—iu — u^) 



(40) 



(41) 



where p is the correlation between W^, W"^. Our general method still applies, however 
the dependence on vq leads to more complicated integrals for R2 which are deserving 
of further study. 



5 Conclusions 

The general picture is that a finite parameter pricing model can be augmented by an as- 
sumption that parameters vary stochastically (hopefully slowly) in time, independently 
of other sources of randomness. Furthermore, one need not assume idiosyncratic errors 
are zero: these too can be modeled by a simple independent random process. These 
assumptions can be built naturally into hedging strategies which minimize risk over a 
given time interval between trades. This risk naturally decomposes into a sum of three 
terms: a term involving changes in the underlier, a term involving changes in the pa- 
rameters and a term involving the idiosyncratic errors. Expectations involving changes 
in the underlier are with respect to the risk-neutral measure, and thus incorporate 
a non-zero price of risk, whereas the remaining probabilities are with respect to the 
physical (historical) probabilities. Decision making for the resulting hedging strategy 
at a given trading time requires knowing how to calculate model prices and generalized 
greeks for fixed values of the parameters (using the analytical formulas, Monte Carlo 
or some other method) and knowing the following data: 

1. the current prices of all exchange traded securities; 

2. historically estimated unconditional variances for the changes in parameters and 
idiosyncratic pricing errors. 

The strategy which results has a number of practical advantages: 
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• It provides a systematic treatment of model errors, hence leading to smaller ab- 
solute hedging errors; 

• Since decision making depends only on current market prices and benchmark 
statistics which arc relatively robust and stable in time, we expect hedging errors 
over successive trading intervals to exhibit only small statistical dependence. This 
means that the probability distribution of the hedging error over longer terms can 
be estimated by the central limit theorem; 

• The method can be implemented efficiently in a wide variety of models, both 
simple and sophisticated, which include pure diffusions, jump diffusions and pure 
jump models. 

It should be observed once more that we are introducing new parameters (here 
V, ©) which describe idiosyncratic errors and the changes in model parameters. Of 
course these themselves are best estimated dynamically in time. Thus the finicky 
practitioner might be lead to introduce further parameters which describe the changes 
in the parameters which describe the changes in model parameters, etc! The logic of 
the paper still applies. 

Further work is suggested along a number of lines: 

• A comparison of the performance of our method relative to the standard discrete 
time hedging strategies; 

• The statistical dependence of successive hedging errors should be studied and 
compared to that of other strategies; 

• A study of whether the strategy is consistent with observed trading patterns in 
the market or whether on the contrary traders use different criteria in decision 
making; 

• In the current paper we avoid discussing transaction costs by assuming a fixed 
schedule of single trades. The work should be extended to include transaction 
costs, and which will allow questions of when it is optimal to trade, and when 
double or higher order trades arc preferable over single trades; 

• The method should also be applied to models for other markets, such as the 
commodity, bond and foreign exchange markets. 

To conclude, the current paper is not intended to present a mathematically con- 
sistent model. On the contrary, it addresses the question of how any mathematically 
consistent pricing model can be extended for use by a rational and prudent practitioner 
who cares about the inevitable errors made by the model. 
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